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ABSTRACT 

We observed 51 Peg, the first detected planet-bearing star, in a 55 ks XMM-Newton pointing and in 5 ks pointings each with Chandra 
HRC-I and ACIS-S. The star has a very low count rate in the XMM observation, but is clearly visible in the Chandra images due to the 
detectors' different sensitivity at low X-ray energies. This allows a temperature estimate for 5 1 Peg's corona of T< 1MK; the detected 
ACIS-S photons can be plausibly explained by emission lines of a very cool plasma near 200 eV. The constantly low X-ray surface 
flux and the flat-activity profile seen in optical Can data suggest that 51 Peg is a Maunder minimum star; an activity enhancement due 
to a Hot Jupiter, as proposed by recent studies, seems to be absent. The star's X-ray fluxes in different instruments are consistent with 
the exception of the HRC Imager, which might have a larger effective area below 200 eV than given in the calibration. 
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1. Introduction 

The star 51 Peg (GJ882 HD 217014) shot to fame in 1995 when 
iMavor & Ouelozl (1 19951) detected an exoplanet in its orbit, the 
planetary parameters being quite unexpected at that time, be- 
cause 5 1 Peg b is a giant planet, located at only 0.05 AU dis- 
tance. The star itself is a G5V star 15.4 pc away from the Sun. 
Its properties are quite similar to the Sun's, since 51 Peg is about 
4 Gyr old and its mass, radius and effective t emperature are com - 
parable to solar values with R = 1.27Rp) (Baines et al. 2008), 
M = 1.1 IMq, T eff ~ 5790K (iFuhrmann et al.lll997l) . However. 
5 1 Peg is a metal-rich star, for which the metallicities given in the 
literature vary o ver a wide range of +0.05 < [Fe/H] < +0.24, 
see for example IValenti & Fischer! d2005l) . Enhanced metallici- 
ties are a common featur e of stars with giant planets dGonzalezI 
I1997L ISantos et al.ll2001l) . 

The activity profile of 51 Peg turned out to be uns pectacu- 
lar. In the Mount Wilson program (Baliunas et al.ll 19951) . which 
monitors the Can H and K line fluxes of main sequence stars, 
the star shows a very low and nearly flat chromospheric activity 
level from 1 977 until 1 989 and a s light drop in 199 and 1 99 1 . In 
the Lowell Observatory program dHall et alj|2007l) . it also shows 
low activity and little variability in Ca fluxes since the beginning 
of the program in 1994. The star was also observed in a 12.5 ks 
ROSAT PSPC pointing in 1992 and detected as a weak X-ray 
source. 

The coronal activity of 5 1 Peg is of interest not only because 
the star is similar to the S un, but also with regard to recent stud- 
ies dKashvap et al.l r2008). which claim stars with close-in giant 
planets to be more X-ray active than stars with far-out ones. 



2. Observations and data analysis 

We observed 51 Peg on two occasions in 2008. A 55 ks XMM- 
Newton was carried out on June 1, 2008, and with Chandra, 
we observed 51 Peg for 5 ks each using HRC-I and ACIS-S on 



December 6, 2008 immediately after each other. The specific ob- 
servation details are listed in Table Q] 



2.1. XMM-Newton data analysis 

The XMM-Newton data were reduced using the Science Analysis 
System (SAS) version 8.0.0. Standard selection criteria were ap- 
plied for filtering the data. In the full-time image obtained with 
the PN detector, the automatic source detection procedure finds 
a faint X-ray source with 32 excess counts at 51 Peg's nominal 
position when using the 0.2-1 keV energy band. This choice is 
motivated by 51 Peg being detected in the 1992 ROSAT PSPC 
pointing as a very soft X-ray source. Because of the weak sig- 
nal, we merged both MOS detectors. In the RGS, no relevant sig- 
nal was present. The PN observation is affected by proton con- 
tamination, therefore we used only time intervals (GTI) where 
the high-energy background averaged over the detector is below 
0.8 cts/s, leading to a PN GTI of 29 ks. 

Since spectral fitting results are not very reliable with this 
low number of counts, we conducted a study in different energy 
bands instead and investigated the recorded counts within the 
source region, a radius of 15" around 51 Peg's nominal position 
for the PN and MOS instruments. The source region size of 15" 
radius was chosen because of the rather broad point spread func- 
tion, which contains 72% (68%) of the photons from a point-like 
source in the PN (MOS) detector. Background counts were ex- 
tracted from source-free nearby regions, which are located on 
the same chip for the MOS detector; for the PN detector, two 
background regions were chosen, one on the same chip as the 
source and one on a neighboring chip. Since 5 1 Peg proved to be 
a very soft X-ray source in the previous ROSAT observation, we 
expect most X-ray photons to be produced from the O vn triplet 
or lines with even lower energies, such as N vn/vi and C vi/v. 
We therefore specified three energy bands for our analysis, con- 
centrating on a band around the Ovu triplet (« 570 eV); the 
detected photons are given in Table [2] Since the EPIC detectors 
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Table 1. XMM and Chandra Observation log of 51 Peg 



Instrument Configuration ObsID Obs. time 

XMM MOS 1 full frame / thick filter 0551020901 

XMM MOS2 full frame / thick filter 055 1020901 

XMM PN full frame / thick filter 0551020901 

Chandra ACIS-S back-illuminated 10825 

Chandra HRC-I imaging 10826 



Table 2. Photons of 51 Peg in XMM and Chandra; see text for 
details. 



Energy PN MOS 1+2 ACIS-S HRC 

(keV) s bs bsbs b 



0.15-0.2 










i+14 
i -0.4 





0.2-0.45 


! 7+ 5.o 

1 '-3.0 


9.1 


7+3.6 
-1.6 


4.9 


6 +3.6 
-1.6 





0.45-0.65 


q+3.9 
°-1.9 


3.5 


6 +3.5 
u -l,5 


3.2 


+2 





0.65-2.0 


j n+4.3 
1 -1,5 


10.4 


21 +72 


23.2 


+2 





0.15-2.0 












... 21+$ 0.6 



have energy resolutions of FWHM^IOO eV, we adopted this as 
minimum bandwidth. 

In the source region we then count the number of photons 
recorded in the various energy bands and detectors, listed in col- 
umn "s" in Table |2] Source counts are giv en with 1 cr confid ence 
limits for low count numb ers according to lKraft et al.l(ll99ll) : for 
a detailed discussion, see lAvresI (12004). The scaled background 
counts (denoted "b") were taken from areas much larger than the 
source region, and thus the error on the background is dominated 
by statistical fluctuations. 

2.2. Chandra data analysis 

For d ata reduction of the C handra observations we used CIAO 
4.1 (iFruscione et al.l 12006]) and applied standard selection cri- 
teria. The analysis of the data was performed in the 0.15- 
1 keV energy band since the back-illuminated ACIS-S chip has 
nonzero effective area at X-ray energies below 300 eV. For the 
HRC imager no energy cuts were used since its energy resolution 
is low. 5 1 Peg is clearly detected in both instruments. 

In the ACIS observation we detect eight photons in the 
source region, a circle with 1.5" radius around 51 Peg's nom- 
inal position. This radius was chosen to contain 95% of the soft 
(< 1 keV) photons from a point-like source. From nearby source- 
free regions in the 150-650 eV energy band we expect only 
0.03 background counts for this area, therefore we attribute all of 
the recorded counts to 5 1 Peg. The spectral resolution of ACIS-S 
is similar to the one of the EPIC detectors (« 100 eV). 

In the HRC-I pointing 2 1 photons were detected in the source 
region over a background of 0.6 photons scaled to the same area. 
At any rate, also the HRC clearly detects 51 Peg. 

3. Results 

3.1. XMM-Newton PN and MOS 

51 Peg shows a photon excess in the 0.2-0.45 keV and the 
O vii band (0.45 - 0.65 keV) in PN and a very weak excess in 
the same bands in the merged MOS detectors. The MOS and PN 
lightcurves show no obvious variability over the whole 55 ks. As 
shown in FigureQ] most of PN's excess source photons have en- 
ergies around 300 eV; another emission feature is present around 



GTI (s) 
55000 
55100 
29000 
4980 
4924 

















IT 


-,-F 


r 




n , , 



100 200 300 400 500 600 700 800 100 200 300 400 500 600 700 800 
PN: Energy (eV) ACIS-S: Energy (eV) 



Fig. 1. Upper panel: effective areas of XMM PN and Chandra 
ACIS-S at energies below 800 eV. Lower panel: detect cell 
counts (solid histogram) of 5 1 Peg in PN and ACIS-S over the re- 
spective backgrounds (dashed; ACIS-S has practically no back- 
ground). 

570 eV, the energy of the O vn triplet. Because of XMM's mod- 
erate intrinsic energy resolution the nominal energies of the de- 
tected source photons cannot be regarded as exact values. From 
the absence of emission features at O vm energies (« 650 eV) 
we can conclude that 51 Peg's corona has an average plasma 
temperature well below 3 MK. 

3.2. Chandra ACIS-S and HRC-I 

All the recorded counts have energies between 150 and 450 eV 
and are distributed quite evenly over the observation time, sup- 
porting a soft, basically constant X-ray source. Let us now in- 
spect the energies of the ACIS-S photons in detail; the CIAO 
software assigns a nominal energy to each recorded photon (see 
Figure [TJ. The eight source photons have energies of 170, 206, 
211, 212, 256, 227, 291 and 428 eV; they are hence very soft 
and obviously none of these photons can be attributed to O vn or 
even O vm emission. This supports our hypothesis of a very low 
plasma temperature evoked by the XMM data. 

The ACIS-S detector is prone to optical contamination, so 
we have to check whether the extremely soft events could be 
induced by optical photons. The threshold for optical contam- 
ination in the ACIS-S detector is at V ~ 7.8 for stars with 
an effective temperature between 5000 and 6500 K; a star this 
bright would cause a bias level shift of one Analog-to-Digital- 
Unit (ADU) of 3.4 eV during the standard 3.2 s time frame for 
the central pixel of the source. 51 Peg's visual magnitude is 5.5, 
so we expect ca. 8 ADUs per time frame. Since the event thresh- 
old lies at 20 ADUs, optical contamination can be ruled out as 
explanation for the detected events. 
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Table 3. X-ray fluxes of 51 Peg with lcr errors, calculated with 
WebPIMMS using a 1 MK thermal plasma model. 



Instrument Flux (0.1-1.0 keV) log L x (0. 1 - 1 .0 keV) 

(erg s citT 2 ) (erg s~') 



XMM PN 


1 2 +LU 


X 


io- 


14 


26.3- 


-26.8 


XMM MOS 1+2 


1 1 +M 
1 - 1 -0.6 


X 


io- 


14 


26.1- 


-26.8 


Chandra ACIS-S 


1 7 +0 - 9 


X 


io- 


14 


26.5- 


-26.8 


Chandra HRC-I 


4 2 +l1 


X 


io- 


14 


27.0- 


-27.2 


ROSAT PSPC 




X 


io- 


14 


26.7- 


-26.8 



Also in the HRC the recorded events are distributed evenly 
throughout the observation time. The intrinsic energy resolution 
of the HRC detector is low so that little information on the spec- 
tral energy distribution can be derived; because the HRC-I ob- 
servation was carried out immediately after the ACIS-S obser- 
vation, we assume that the 21 detected HRC-I source photons 
have similar energies as the photons in the ACIS-S detector. 

3.3. Determination of coronal temperature 

To estimate the coronal temperature, we evaluate the 
temperature-dependent hardness ratios of several energy bands, 
viz. HRpn = Hppf/Sptf with Hp^ and S pn covering 0.45- 
0.65 keV and 0.2-0.45 keV for PN; for ACIS-S, we use 
HRacisi - H AC is i/S acis i and HRacis 2 - Hacisi/Sacisi 
with the energy bands Hacis i- 0.25-0.45 keV, Sacis i- 0.15- 
0.25 keV, Hacis 2- 0.45-0.65 keV and Sacis 2- 0.25-0.45 keV. 
The energy bands for ACIS-S were chosen to obtain quite an 
even distribution of source photons. 

We derive the temperature-dependence of the hardness ra- 
tios with Xspec vl2, using the instrument responses and ef- 
fective areas as shown in Figure Q] and simulating spectra for 
a one-temperature thermal plasma model with solar abundance 
in the temperature range between logT = 5.6 and 6.6 (see 
Figure 13. The possible plasma temperature of 51 Peg is then 
constrained by the observed hardness ratios for ACIS-S and 
PN. Assuming lcr errors for low count numbers as defined in 
Table El we find that HR PN = 0.6^, HR AC is 1 = 0.6^ ° and 
HRacis 2 = 0^. The observed PN ratio yields the tempera- 
ture limits 5.85 < logT < 6.3. The ACIS-S ratios yield a 
lower temperature limit from HR2 and an upper limit from HRy. 
5.8 < log T < 6.05, so that the likely temperature range for 
51 Peg's corona is 5.85 < log T < 6.05. 

Since there are virtually no background photons in ACIS-S, 
we can use the energies of the eight recorded source photons 
to validate our temperature constraints by identifying the most 
likely emission lines of their origin. Strong emission at temper- 
atures near 1 MK comes from the Si and S emission line com- 
plexes around 200 eV, some strong Si lines around « 230 eV, 
the C v triplet around 300 eV and the N vi triplet around 426 eV. 
These emission lines match well with the observed photons, 
which can be considered as a rough plausibility check. So, both 
the temperature constraints from hardness ratios and the identi- 
fication of possible emission line complexes point to a plasma 
temperature of T <, 1 MK. 

3.4. Comparison with ROSAT data 

To investigate the long-term X-ray behavior of 51 Peg, we re- 
analyzed a 12.5 ks ROSAT PSPC observation from December 




og T log T 

(c) HR PN (d) HR PSPC 

Fig. 2. Hardness ratios as a function of temperature (see text for 
details). Observed ratios are given as crosses with lcr error. 



1992. Using the conversion factor ecf = (5.30 x HR { + 8.31) x 
10 2 erg count 1 cirT 2 with HR\ being the hardness rati o HR\ = 
(H - S)/(H + S) (S: 0.1-0.4keV, H: 0.5-2.0 keV) (ISchmitj 
1 1997b . the luminosity derived from the observed count rate of 
7 cts/ks is \ogLx = 26.75 erg s _1 . The observation is split in two 
parts, interrupted for ca. 70 ks. 51 Peg emits mainly soft X-ray 
photons below 500 eV (HR[ as -1.0) and showed a low and 
constant X-ray activity level without any obvious variability. 

We can constrain the coronal temperature of 5 1 Peg in this 
observation in the same way as for the XMM and Chandra 
pointings. We use the temperature-dependent hardness ratio 
HRpspc = H p S pc /S pspc with the energy bands H P sp C : 
0.1-0.3 keV, Spspc- 0.3-0.65 keV. The observation yields 
HRpspc = 0-12*oo6 fr° m which lower and upper limits for the 
temperature can be derived, namely 5.65 < logT < 6.05 (see 
Figure|2]l. This again leads to a temperatur estimate of T < 1MK. 

3.5. Consistency of measured count rates and fluxes 

The observed count rates in XMM PN and MOS and Chandra 
ACIS-S are, considering the low photon numbers and therefore 
large statistical fluctuation, in reasonably good agreement. That 
ACIS-S detects no photons in the O vn band is no surprise given 
the shorter exposure in ACIS-S and its smaller effective area in 
that energy range (see Figure [TJ. At energies below 300 eV, we 
find that the PN and MOS counts numbers are smaller than ex- 
pected from what we see in ACIS-S. This might be explained 
by statistical fluctuations, errors in the effective area determina- 
tion or energy redistribution effects in the CCD detectors (EPIC 
"low-energy shoulder"). Given these uncertainties for very low 
energies, we use only the O vn counts (0.45-0.65 keV) of PN and 
MOS for our flux calculations and then extrapolate the flux to a 
common energy range of 0.1-1 keV for comparison. For the 
other instruments, we use 0.65 keV as the upper bound of the 
energy range and their low-energy sensitivity limits as the lower 
bound (0.15 keV for ACIS-S/HRC, 0.1 keV for PSPC) and then 
extrapolate to the common energy range. 

The fluxes normalized to the 0.1-1 keV energy band and the 
corresponding X-ray luminosities are consistent within lcr errors 
except for the HRC-I flux, which seems to be larger. The count 
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mean S index from 1977-1991 
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1985 1990 1995 2000 2005 2010 

Year 

Fig. 3. X-ray luminosity (crosses) and S index seasonal mean of 
51 Peg's Can H and K line flux; reliable Ca data given as filled 
triangles with typical standard deviations of < ±0.005, open tri- 
angles are less reliable values derived from only a small number 
of observations. 

rate measured by the HRC instrument is higher by a factor of ca. 
2.5 compared to the ACIS-S count rate. The nominal effective ar- 
eas of the two instruments are very similar at low energies, with 
the HRC having somewhat larger effective area below 200 eV 
(AA * 10 cm 2 or 17% at 200 eV, 44% at 150 eV). The additional 
counts might arise from photons at these energies, but consider- 
ing the small difference in effective areas, it does not seem likely 
that this is the case for all excess HRC photons. This mismatch 
is further validated by comparing ACIS-S and HRC count rates 
with WebPIMMS: assuming a thermal plasma with solar abun- 
dances and T = 1 MK, 8 counts in the 0.15-0.65 keV energy 
band in ACIS-S translate into 9 expected counts in the same en- 
ergy band in HRC-I, which is obviously inconsistent with the 
21 recorded HRC photons only 15 minutes after the ACIS-S ob- 
servation. The photon count estimate changes by < 20% if one 
assumes a plasma temperature of 0.8 or 1.25 MK, so a slightly 
different plasma temperature does not cure the substantial mis- 
match in the count rates. 

Mismatches between HR C and XMM count rat e s have been 
repor ted before for a Cen (Ro brade et alj 12005; Avre s et all 
2008). This mismatch between almost simultaneous HRC and 
ACIS-S count rates can be explained reasonably by two pos- 
sibilities: either the effective area of the HRC at low energies 
is underestimated in the current calibration or the effective ar- 
eas of XMM MOS and PN (while using the thick filter) as well 
as Chandra ACIS-S are overestimated in that energy range. A 
detailed cross-calibration effort, preferably with a soft coronal 
source, could help to resolve any systematic errors in the effec- 
tive areas of the instruments. 

4. Discussion 

4.1. Low activity - 51 Peg a Maunder minimum candidate? 

We found 5 1 Peg to be a rather constant, weak and soft X-ray 
source over the last 15 years. Another available long-term ac- 
tivity indicator is Can. In the Can H and K line flux moni- 
toring program s carried out at the Mount Wilson an d Lowell 
Observatories dBaliunas et all 119951; lHall et alj 120071) . 51 Peg 
was found to have a very low chromospheric activity level 



((R HK ) - -5.01, S mw ~ 0.16). In Figure[3]we plot the star's 
Mount Wilson S index measured since 1996 together with the 
average of older data. Clearly, the overall chromospheric activ- 
ity is low, with some variations in the older set of data which is 
also seen in the more recent observations. Apart from one data 
point which is derived from a very small number of observa- 
tions, 51 Peg's S indices are at the lower end of or even below 
the Sun's respective dat a during a solar minimum (data taken 
from Bali unas et alj (Il995l) ). Other stellar properties like radius, 
mass, age and effective temperature are similar to the Sun's re- 
spective parameters. 

The steady low-activity behavior of 51 Peg's Can H and K 
line fluxes is also reflected by its X-ray properties. Compared to 
estimates for the solar X-ray luminosity in the ROSAT RASS 
band (0.1-2.4 keV) during a solar cycle (lJudge et alj 120031) . 
51 Peg's luminosity is also at the lower end of the Sun's val- 
ues. The ratio of the star's X-ray to bolometric luminosity is 
also rather low with Lx/Lboi = 1 X 10~ 7 . The X-ray sur- 
face flux of F to M stars was shown to be constrained at the 
lower end by the surface flux level of a solar coronal hole; 
^x(hole) ~ 10 4 el "g s 1 cmT 2 for the ROSAT and Chandra en- 
ergy band, which tr anslates to » 10 38 erg s _1 cmT 2 for XMM' s 
0.2-12 keV band dSchmittll 19971) . 51 Peg's surface flux, calcu- 
lated from the ACIS-S data, is one of the lowest so far detected 
with log fx = 10 3 7 erg s _1 cmT 2 ; the coronal hole surface flux 
seems to be a good description of this star's X-ray flux, with 
regards to the flux level as well as the plasma temperature. 

There has been some discussion on how to identify a 
Maunder minimum (MM) star over the last years. The orig- 
inal criterion o f chromospheric ac tivity levels (R HK ) - -5.1 
was derived by iHenrv et al.l (1996), but relied on a stellar sam- 
ple contaminated with evolved stars, which have significantly 
lower chromospheric activity levels compared to main sequence 
stars. I Wright! d2004l) reanalyzed these data, excluding evolved 
stars with luminosities more than one magnitude larger than the 
Hipparcos average main sequence for the respective B — V value, 
and found that most stars previously identified as MM candi- 
dates are evolved stars and therefor e not comparab l e to th e Sun's 
Maunder minimum state. This led Judge & Saar] d2007l) to the 
question if the minimum (R HK ) level for main sequence stars to 
qualify as an MM candidate should be higher than -5 . 1 , and also 
to consider flat-activity time profiles and U V- and X-ray data 
to identify MM candidates. A recent study by lHall et all (G009) 
suggests that minimum levels of R HK depend on stellar metal- 
licity, with metal-poor stars from the examined sample having 
a higher minimal R HK - In this picture, 51 Peg as a metal-rich 
star still has low chromospheric activity as measured by R HK , 
but this alone does not necessarily qualify it to be a Maunder 
minimum c andidate. However, as recent results show ( Hal l et ail 
120071 12009). the absolute magnetic excess flux AJ^hk seems to 
be a more reliable indicator for stellar activity than R HK - In terms 
of this quantity, 51 Peg's activity level is even lower compared 
to the quiescent Sun than indicated by R HK or the S index, sup- 
porting our interpretation of 5 1 Peg as being extremely inactive. 

The strongest line of evidence for 5 1 Peg being a Maunder 
minimum candidate is its flat activity p rofile as seen over 
decades in the Mount Wilson program (Baliunas e t alj 1 1 995h 
and in observations at Lowell Observatory (lHall et ai]|2007l) . as 
well as the extremely low X-ray surface fluxes, which have not 
changed significantly since the 1992 ROSAT obs ervations. That 
51 Peg is a slow rotator with a 30 - 40 d (Bal iunas et al.l 
ll996l;lMavor & OueloJ 19951) fits the picture, making 5 1 Peg the 
first MM candidate star with a close-in giant planet. 
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A statistical analysis of the X-ray luminosi ties of planet- 
bearin g host stars has recently been conducted (Kashyap et al. 
2008). Its authors claim that stars with close-in giant planet, such 
as 5 1 Peg, are on average X-ray brighter by a factor of two com- 
pared to stars with far away planets. Apparently, 51 Peg's over- 
all activity is not enhanced by the presence of its Hot Jupiter. 
However, at a distance of order of 50 R j„ p only a weak interac- 
tion between an inactive star and its planet might be expected. 



Santos, N. C, Israelian, G., & Mayor, M. 2001, A&A, 373, 1019 
Schmitt, J. H. M. M. 1997, A&A, 318, 215 
Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141 
Wright, J. T. 2004, AJ, 128, 1273 



5. Conclusions 

We have detected X-ray emission from 5 1 Peg in a 55 ks obser- 
vation with XMM-Newton and 5 ks observations with Chandra 
ACIS-S and HRC-I each. The detection of 51 Peg with a low 
count rate in the XMM pointing and the clear source signal in 
the much shorter Chandra observations can be explained by the 
different effective response of the detectors at low energies and 
5 1 Peg having an extremely cool corona. Our main results are 
summarized as follows: 

1 . 51 Peg shows weak emission in the O vn triplet and emission 
around 200 eV which can be explained most likely by cool 
silicon emission lines. 

2. A coronal temperature of < 1 MK is consistent with the de- 
tected hardness ratios in different energy bands in both the 
XMM and the Chandra pointing as well as in the ROSAT ob- 
servation carried out 16 years earlier. 

3. The Chandra HRC-I count rate is higher than can be ex- 
plained by differences in the effective areas of HRC and 
ACIS-S; HRC's effective area might be larger at low ener- 
gies than given in the calibration so far. 

4. The constant and very low surface X-ray flux level together 
with a flat-activity behavior in chromospheric Can H and 
K line fluxes suggests 51 Peg to be a Hot Jupiter-bearing 
Maunder minimum candidate. 
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